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Abstract 

In this paper, we partially answer open questions about the convergence of overlapping 
Schwarz methods. We prove that overlapping Schwarz methods with Dirichlet transmis- 
sion conditions for semilinear elliptic and parabolic equations always converge. While 
overlapping Schwarz methods with Robin transmission conditions only converge for semi- 
linear parabolic equations, but not for semilinear elliptic ones. We then provide some 
conditions so that overlapping Schwarz methods with Robin transmission conditions con- 
verge for semilinear elliptic equations. Our new techniques can also be potentially applied 
to others kinds of partial differential equations. 

Keywords: Domain decomposition, Schwarz methods, semilinear parabolic equations, 
semilinear elliptic equations. 
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1. Introduction 

The Schwarz domain decomposition methods are procedures to solve partial dif- 
ferential equations in parallel, where each iteration involves the solutions of the original 
equation on smaller subdomains. The alternate method was originally proposed by H. 
A. Schwarz [24] in 1870 as a technique to prove the existence of a solution to the Laplace 
equation on a domain which is a combination of a rectangle and a circle. The idea was 
then used and extended by P. L. Lions [lf|, [ltjj . [l7j to parallel algorithms for solving 
partial differential equations. Since then, many kind of domain decomposition methods 
have been developed, to improve the performance of the classical domain decomposition 
method. However, the convergence for domain decomposition methods still remains an 
open question. 

Many techniques have been developed to prove the convergence of classical Schwarz 
methods, or Schwarz methods with Dirichlet transmission conditions. One of the first 



techniques, used by P. L. Lions in 15j], is the iterated projections for linear Laplace 
equation and linear Stoke equation. The idea is to prove that classical Schwarz methods 
for these equations are equivalent to sequences of projections in Hilbert Spaces. In the 
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same paper, P. L. Lions also showed that the Schwarz sequences for nonlinear monotone 
elliptic equations are related to classical minimization methods over product spaces and 
proposed to use Schwarz methods for evolution equations. This idea was then used by L. 
Badea in [l[ to prove the convergence of classical Schwarz methods for nonlinear mono- 
tone elliptic problems. 

Following the pinoneering work of P. L. Lions, in the papers Q, Q, E. Giladi, 
H. B. Keller, A. Stuart and M. Gander used Fourier and Laplace transforms, together 
with some explicit calculation to study classical Schwaz methods for some 1-dimensional 
evolution equations, with constant coefficients. Later, by using a maximum principle 
argument, M. Gander and H. Zhao proved that classical Schwarz method converges for 
the n-dimensional linear heat equation [8]. 

Another technique to study the convergence of classical Schwarz methods is to use the 
idea of upper-lower solutions methods, with initial guess to be upper or lower solutions 
of the equations. This special class of domain decomposition methods with monotone 
iterations has been studied by S. H. Lui in (ljj) . (20| . 21 1. Although many techniques 



have been developed to study the convergence problem of classical Schwarz methods, the 
problem with nonlinear equations in n-dimension and general multi-subdomains is still 
open. 

A new class of Schwarz algorithms, in which Dirichlet transmission condition is re- 
placed by Robin ones, has been studied recently in order to improve the performance 
of classical methods. The new algorithms are called optimized Schwarz methods since 
there are some parameters we can optimize to get faster algorithms. In 1989, P.L. Lions 
(see [lfjj], [l?]]) established the convergence of nonoverlapping optimized Schwarz methods 
with Robin transmission conditions by using an energy argument. Later, J. D. Benamou 
and B. Depres in Q used this technique to study the convergence of nonoverlapping op- 
timized Schwarz methods for Helmholtz equation. Energy estimates have then become a 
very powerful technique to prove the convergence of nonoverlapping optimized Schwarz 
methods with Robin transmission conditions (see (HI). 

However, the convergence problem of overlapping optimized Schwarz methods, even for 
linear problems, still remains an open problem up to now. J.-H. Kimn [131 ] . proved 
the convergence of an overlapping optimized Schwarz method for Poisson equation with 
Robin boundary data, 

f -Ait = f in.il, 

\ $%+jm = gondCl. 

He proved that there is an pq > such that the Schwarz iterations with Robin transmis- 
sions conditions converge for any Robin parameter < p < pa . In 0, S. Loisel and D. 
B. Szyld extended the technique of J.-H. Kimn for the following equation 

( -V(aVu) + cu = f in Q, 

[ u = on <90, 

where a is a C 1 -function and c is positive and belongs to L°°(Cl). The same constant 
p is kept for all transmission operators and some conditions on the boundaries of the 
subdomains are then imposed. 

A proof of convergence based on semi-classical analysis for overlapping optimized Schwarz 
methods with rectangle subdomains, linear advection diffusion equations on the half plane 
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was given in p3| . 



Another technique is to use Fourier transform. This technique cannot be used to study 
the convergence of Schwarz methods for nonlinear problems and for general subdomains, 
but convergence rates can be obtained. Changing the boundary conditions will change 
the values of the convergence rates and then improve the performance of the algorithms, 
which proposes a new problem: the problem of optimizing the convergence rates. In [l2j ]. 

, [f| the authors showed that the problem of optimizing the convergence rates is in fact 
a new class of best approximation problems and suggested a new method to solve it. 

In this paper, we present convergence proofs of overlapping classical and optimized 
Schwarz methods for elliptic and parabolic semilinear equations, in general forms, for 
general multi-subdomains. We prove that Schwarz methods with Dirichlet and Robin 
transmission conditions always converge for parabolic equations; since with parabolic 
equations the time variable can be controlled easily. However, Schwarz methods with 
Robin transmission conditions do not converge for elliptic equations, while classical 
Schwarz methods always converge. We can see from Remark 3.1 that given a Schwarz 
algorithm with a specific Robin transmission condition, there exists a class of elliptic 
equations where the algorithm is unstable. A condition of convergence is then supplied: 
Schwarz methods with Robin transmission conditions for elliptic equations will converge 
if we multiply Robin parameters by a number large enough, and this can also be seen 
from Example 3.1. The techniques used in our proofs can also be used to prove the 
convergence of Schwarz methods for many other kinds of partial differential equations. 

The paper is organized as follows. 

• Section 2 is devoted to the convergence properties of Schwarz methods for semilin- 
ear parabolic equations. Section 2.1 gives the definition of the Schwarz algorithms 
for semilinear parabolic equations, and states the two theorems of convergence. 
Theorem 2.1 announces that classical Schwarz algorithms always converge with 
semilinear parabolic equations and its proof can be found in section 2.2. Theo- 
rem 2.2 is about the convergence of Schwarz algorithms with Robin transmission 
conditions and the proof is then given later in section 2.3. 

• In section 3, we discuss the convergence properties of Schwarz methods for semilin- 
ear elliptic equations. Definitions of the algorithms, the two convergence theorems 
3.1, 3.2, together with the counterexample 3.1 is announced in section 3.1. Section 
3.2 and 3.3 contain the proofs of the two theorems. 

2. Convergence for Semilinear Parabolic Equations 

We introduce the abbreviation = dx d dx , dt = and di = and consider 
a general semilinear parabolic equation 



I. u(x, 0) = g{x, 0) on f2, 

where Q is a bounded and smooth domain in R n . The coefficients dij, bi, c are functions 
of the space variable x, with the following properties 




dj(dijdiu) + X)i=i bidiii + cu = F(x, i, u) in x (0, oo), 



(2.1) 
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(Al) The functions a itj , h, c are in C 2 (R"). 

(A2) For all i,j in {1, . . . , /}, a^j(x) = aj^(x). There exist strictly positive numbers A, 
A such that A — (a,ij(x)) > XI in the sense of symmetric positive definite matrices and 
\ciij(x)\ < A in Q. 

(A3) g is in C 2 (R™ +1 ) and and F is uniformly Lipschitz in the third variable, i.e. there 
exists C > 0, such that 

V t e K, V x e K T \ \F(x,t,z) - F(x,t,z')\ < C\z - z'\, V z, z' e K. 
With Conditions (Al), (A2) and (A3), Equation (|2.ip has a unique bounded solution it 
in C 2 ' 1 ^ x (0, oo)), i.e dijd t u belongs to C(Q x (0, oo)) for all Lj in {1, . . . , n}. The 
proof of this result can be found in some classical books like 0] , [3] ■ 

The domain is divided into I smooth overlapping subdomains {f2; such 
that 

u? =1 n, = n ; 

{dSli\dSl)n(dSl v \dSl) = 0, V 1,1' G {1,...,/}, l + l'; 

and 

viG{i,...,/},vr,i"G J h i"^i', n,/nn,// = 0, 

where 

j/ = n Qi ^ 0}. 

For any I in J, for I' E Ji, Fiji is the set (<9f2;\<9f2) n fi;'. 




Figure 1: A good way of dividing Q 



Remark 2.1. Figure 1 gives an example which satisfies our assumptions about the way 
we divide Q into several subdomains. In Figure 2, since there is an overlapping area 
between the three subdomains, this way of dividing ft does not satisfy our conditions. 
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Figure 2: A bad way of dividing Q 



The Schwarz waveform relaxation algorithm solves / equations in / subdomains in- 
stead of solving directly the main problem (|2.1[) . The iterate #fc in the Z-th domain, 
denoted by uf, is defined by 



where the transmission operator *Bij> is either of the Dirichlet type or of the Robin type. 
Each iterate inherits the boundary conditions and the initial values of u 



A bounded initial guess it in C°°(fi x (0, oo)) is provided, i.e. at step 1 Equations (|2.2[) 
are solved 

KiMul = u° on T hV x (0,oo),W G J/. 
We assume also the compatibility condition on u 

Bi, r g(.,0)=u°(.,0) onr V) Vl' e J h 

Denote by ef the difference between itf and u, and subtract Equation (|2.2p with the 
main equation (|2.1[) to obtain the following equations on ef 




+ cuf = F(t, x, uf), in x (0, oo), 



(2.2) 



uf = g on (dfli n dfl) x (0, oo), uf (., 0) = 0) in 



f ftef - £™ J=1 djiaijdie* ) + ElLi 6 ^t e f + ce f = «?) ~ s, «) m ^ x (0, °°) 




ef (., .) = on (<9f2 ; n 90) x (0, oo), ef (., 0) = in 
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2.1. Classical Schwarz Methods 

Consider the classical Schwarz waveform relaxation algorithm with Dirichlet trans- 
mission conditions 03;^' = Id. By induction, each subproblem (|2.2j) in each iteration has 
a unique solution in C 2 ' 1 ^ x (0, oo)) then in L 2 (0, oo, fl" x (0)) n L°°(n x (0, oo)) also. 
Consider (|2.3I) and let g, / be bounded and strictly positive functions in C°°(R™,R) and 
C°°(R,R). Define 

*t(x,t) := (ef) 2 g(x)f(t). 

Since ef belongs to L 2 (0, oo, H\n)) U L°°(rj x (0,oo)), $f belongs to L 2 (0, oo, H 1 ^)). 
Let c, be 6, + 5Z" =1 ^i.jdjgg -1 , then Cj € L°°(S1; x (0, oo)), and define the following 
operator 

n n 

Lemma 2.1. 7?t, eac/i subdomain for each iterate k 

m £/ie distributional sense, i.e. for all ip in Hq(Q) and <p > a.e. on Q, 

a.e. in (0, oo). 
Proof. Define the operator 

n 

-Cfijo -=d t - E dj(a,idi). 
A lengthy but easy computation then implies that 

(n \ n 

fte? " E ,V - KA.?4') ef 5/ - E Zatjfafaetgf - (2.5) 

n / n 

- J2 ^jdiefefdjgf + (ef ) 2 - ^ KAiS + d iaid d j9 )f + gf 

Thanks to Equation (|2.3[) , and the lipschitzian property of F, the first term in (12.51) can 
be estimated in the distributional sense 

2 ^ef - E ^(^^ef)j ef 5 / (2.6) 

= 2 5 /ef ^F(t, uf ) - F(t, as, u) - E b i d i e l ~ ce ^ 
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re? 



< 2gf(ef) 2 (C+\\c\U-J2^d l efgfe 

i=l 
n 

< 2gf(e*) 2 (C+ \\c\U (ft*? - (ef) 2 ^/) 

i=l 

n / n ^ 

< + (ef) 2 2 3 /(C+|| C || 0O ) + ^6 l 9, g / 
i=i V »=i y 

Since £ido is elliptic, the second term in (I2.5[) is negative. 
Moreover, the third term in (12.50 can be transformed into 



- E 4a i , i S i efef9 i3 / E ^d.gg^d^ + E 2^%^ L /(ef) 5 
»j=i <ii=i *,j=i 

Combine ([2~5]l. and to get 



W*f) + E | 6 ' + E ^ j9 g- 1 ] ft*? < (ef ) 2 97l, (2.7) 

2=1 



where 



£01 = 



/ ftjff ~ % d j gd l g f ^ ftg 

> {-au — djdij \-2dij— +- + 2C+ c oo +) 6* — 

id^i 9 9 9 9 I j-[ 9 



Notice that 9Jt has the form of fg(y- + G(g)). Choosing / such that — y- large enough 
(for example, / = exp(— at), where a is a large positive constant), since the other terms 
are bounded in the bounded domain Qi x (0, oo), then 9Jt < 0. The nonlinear equation 
(|2.3p has been transformed into the following linearized inequation of *? 

n n 

ft*? ~ E <>,i", ,'>M: + E ^ °- ( 2 - 9 ) 

□ 

Theorem 2.1. Consider the Schwarz algorithm with Dirichlet transmission condition, 
suppose that f(t) is a strictly positive and continuous function satisfying — min tg (o,oo) yjry 
is sufficiently large, we get the geometrical convergence 

lim max ||(wf - w) 2 /WIU~(o ! x(o,oo)) = 0. 

Remark 2.2. In the above theorem, if f is chosen to be exp(— at), then when a is large 
enough, 

. /'(*) 
- mm — — 

te(0,oo) f(t) 
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is large enough. In this case, the limit 

lim max ||(uf -M) 2 /WllL~(n i x(o,oo)) = 

implies the almost everywhere convergence of the sequence {uf} to u on f2/ x (0,oo). 

Remark 2.3. In the proof, if aij, hi, c depend both on t and x, the convergence result 
in the theorem remains true. 



Step 1: Construct estimates of the errors ef from Inequation (|2.9p . 



Proof. The proof is divided into two steps 
Define 

M = esssup $f (#,*), (2.10) 
ar2,x[o,oo)uO;x{o} 

we will prove that the maximum principle holds, i.e. M > $f a.e. on fli x (0, oo). Define 
the function 

w = ($f - M)+ = max{$f - M, 0}. 
Since w E L 2 (0, oo, H^(Q t )), then 

n n 

d t w — <),■(!, . ,<),<!■ : + \^ CjSjiu < 0. (2-11) 

i,j=l «=1 

To prove that M > $f a.e. on 51; x (0,oo), it suffices to prove that w = a.e. on 
Hi x (0, oo). 

For < p < oo and < q < oo, and for < t\ < r 2 < oo, define | \h\ \ Ti ,t 2 . p. g = 
II^IU^^^^LPfo,))' for h e L q (T 1 ,T 2 ,L p (fli)). If ti = 0, denote ||/i||o,T 2 , P ,g by \\h\\ T2tPtq . 

Let x{t\,T2) be the characteristic function of the open interval (ti,T2), where < 
n < r 2 < oo and set 99 = u>x- Since w € L 2 (0, T, i?g (£!/)) nL°°(f2 ; x (0, 00)), it is evident 
that <p e L 2 (0, 00, H&(Qi)) nL°°(Qi x (0,oo)). 

Use ip as a test function for (12. lip 



dtwwdxdt + / / aijdiivdjwdxdt + / / 'y^CjdjWwdxdt < 0.(2.12) 

Jti Jf2; jj- = i •'Ti JOi i=1 

Equation (|2.12[) and Conditions (Ai) and (A2) imply that 

f ^-dx + I I \\Vw\ 2 dxdt< I I M^Vwlwdxdt, (2.13) 



where Mi is a positive constant. By the Cauchy inequality, the right hand side of (|2.13[) 
is bounded by 

AI 6 f T2 I AI 

-^-\Vw\ 2 dxdt + / / —±w 2 dxdt, (2.14) 
2 J Tl Jn, 2e 
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where e is a small positive constant. 
For e to be Equality (|2.14l) implies 



f W 2 t ~ T2 r T 2 r \ rT2 r 

/ — dx + / -\Vw\ 2 dxdt < / / M 2 w 2 dxdt, (2.15) 
Jsii 2 t=Tl Jti Jn, 2 Jn JUi 



with M 2 = 

Denote -X'(t) = J w 2 (x,t)dx, and let T2 be i in the interval 3 = (ti,tx + 8), the 
previous estimate infers 

X{t) + X\\Vw\\ 2 Tit22 < M 2 8{supX(t)} + X{ Tl ). (2.16) 

a 

Choosing 8 such that M2<5 = h, the fact that sup 3 X{t) < 2X{t\) then follows. Since the 
inequality is true on any time interval with the length of <5, and X(0) = 0, then X(t) = 
for a.e. t in (0, oo). Hence w — for a.e. t in (0, oo). 
We have just proved that 

(ef(x,t)) 2 g(x)f(t) <max( esssup (ef(x,t)) 2 g(x)f(t) ) , (2.17) 

for all Z in /, for a.e. (x,t) in x (0, oo); for any strictly positive functions g, f in 
C*°°(R",R) and C*°°(M,R). 
Step 2: The convergence of the algorithm. 
Denote 

E k = max ( esssup (ef) 2 f(t) ) . (2.18) 



(z,t)6(fi(X(0,oo)) 



From (|2.17p comes that for every V in Ji and for a.e. (x,t) in Tij> x (0,oo) 



(ef(x,t)) 2 5 ( 2 ;)/(i) < max esssup (efr^x.t)) ff(ar)/(t) , (2.19) 

Vr,,,»x(o,Do) / 



that implies 

(e?(z,*)) 2 /(i)< ^T, m ^ ( eswup (ef,r 2 (x,i)) 2 <?(z)/(i) ] . (2.20) 

Since r^// lies inside J7/< , choose g such that there exists a constant satisfying 

1 > M u > VC G I> and VC G U^^IV,,,,, 

5(0 

and this implies that for all I' in J;, for a.e. (x,t) in T;^' x (0, oo) 

(ef(x,t)) 2 /(t) < Af 3 . z max ( esssup (ef,r 2 (x, i)) 2 /(*) ) < M 3tl E k ~ 2 . (2.21) 
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Choose g to be the function 1, (I2.17|) yields 

(ef(x,t)) 2 f(t) < max ( esssup (ef{x,t)) 2 f(t)\, VZ' £ J;, a.e. on f2;. (2.22) 

\r M , x (°.°°) * / 

The estimates (|2.21[) and (|2.22l) imply the existence of a constant M4 smaller than 1 
and satisfy 

E k < M 4 E k - 2 , (2.23) 
which shows that the errors converge geometrically 

lim E k = 0. 

The theorem is proved. □ 

2.2. Optimized Schwarz Methods 

The optimized Schwarz waveform relaxation algorithms are defined by replacing the 
Dirichlet by Robin transmission operators 

n 

?Bi,l'V= ^2 aijdivnij'j +pi,i>v, 
«.i=i 

where riu'.j is the j-th component of the outward unit normal vector of IY/'; pij/ is 
positive and belongs to L°°(IY//). By induction, each subproblem (|2.2p in each iteration 
has a unique solution in L 2 (0, 00, H 1 ^)) and the algorithm is well-posed. 
Let / be a function in L 2 (0, 00), define 



f(x) exp(—yx)dx. 
Now, define for a fixed positive number a 



a' + l / />oo \ 2 

f(x)exp(-yx)dx) dy 



\f\a = SUp 
a' >a 

and 

L*(0,oo) = {/ : /eL 2 (0,oo),|/| Q <oo}. 

Then (L„(0, 00), |.| Q ) is a normed subspace of L 2 (0, 00). 

Consider Equation (|2.3p . let gi be a function bounded and greater than 1 in C°°(M™, R), 
a be a positive constant, and define 



$f(x) := Qf ef exp(-at)dtj gi{x), 
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then $i(x) belongs to H 1 ^). 

Let B\ and C l be functions in L°°(R n ) 



defined in the following ways: 



j=i V yz y 



2d t gid 3 gi dig d ltJ gi 

(gi) 2 g 91 



Define 



(2.24) 



+ {jT [(f + c) ef - F(u?) + F(u)] exp(-crf)di J , 9i . 



(2.25) 



ft is possible to suppose a to be large such that C l belongs to (j,a). 

Lemma 2.2. Choose gi, gi> such that Vgi = Vgi' = on Fiji and ^ > 1 on ^i,i>, for 
all V in Ji . <3?f is then a solution of the following equation 

j £ m ($f) = 0, infij x (0,oo), 

1 /3 Z ®,,,, ($f ) =»i,i' onI V x (0,oo),Vr g J;. 

w/iere pi = ^ on Fiji , for all I' in Jj. 

Proof. A complicated but easy computation leads to 

n 

/■oo ( n \ 

/ 9 * e i fe _ E d i( a iJ d i e l) exp(-at)dt 

/"CO \ / n 

y exp(-at)dtj I a,. ,[(>,. .w • il<' , ,<Kgn 



(2.26) 



n r /*oo 

^ ffli,j / ef (d l e^djgi + dje^d l gi)eyLv{-at)dt 
i,i=i Uo 



That implies 



(2.27) 
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J ef cxp(-at)dtj I ■ n,,,i),, i gi + djdijdigi) 

y ( ~ E bidi6 i ~ ce ' fe + F (^ fe ) ~ F ( u n e M-at)dt 

ef (d l efd j gi + djefdigi) exp(-at)dt 



91 



-E 



Therefore 
= 



n n n / o o \ 

- ^ »j(aijft*f)+x; w + E «ij (ft^^+flj*?^) 

i,i=l 1=1 i,j=l ^ gi gi ' 

+ (?+E ^(-^^+^^-^^)-E^) *? + 



i=l 



+ {jT [(f + c) ef - F(u?) + F(u)] exp(-ai)^J 5; . 

Now, consider Robin transmission conditions on the boundary ry and notice that 
\7gi = Vgi' — on Fiji , the transmission conditions can be reported on $f 

n 

®U'(*?) = E ^SjWj+PU'*? (2-28) 

„oo / n \ n 

= / E a hi n l,l'J d i e l +Pl,l' e i exp(~at)gi+ ^ a^-ny^-d^ef 
E a h3 n l,l'J d i e l'~ 1 +Pl,l' e i~ 1 cxp(-at)gi 

Vj=i y 

n \ 



□ 



Theorem 2.2. Consider Schwarz algorithms with Robin transmission conditions. There 
exists a constant a such that for a to be greater than a 



lim £ / \ef\ 2 a dx = 0, 

1=1 Jit l 



Proof. For all I in {1,/}, denote by to be the open set f2;\Ui' e j,0/'. Let </?f be 
functions in H 1 ^) and yf +1 be functions in if 1 ^/) for all I in / such that <^f +1 = </?f 
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on Fij/ for all I' in J;. Now, use <^f +1 and tpf as test functions for (|2.25p . and take the 

sum with respect to I in {1,/} the integrals £,m($>i)(pf and £i_R($f +1 )<,sf +1 , then 

I ( „ n n . . 

E W- E ",^*f/<U'/.<- • E /- B l A*i<Pi<b+ I C l $?rfdx 

1=1 \Jci'i,j=l i=i'' n i J n i 

E / Pi'^irfda 
re./, Jt i',i 





r r 00 






/J 










Jo 







(2.29) 



1=1 [• / "IU=1 i=l 



- + cj ef +1 - F(uf +1 ) + F(«)J exp(-crf)cft }> gitf +1 dx 



In the above equality, choose </?f +1 to be then there exists <y9; fc , such that 

<^ = <^|"/ +1 on r^;/ for all /' in J;; moreover, 

llvfllffitfi,) < ^ E ll^?' +1 IUi(fV) and llvfll^cno < C E H^IU'PV)' 
where C is a positive constant. 

With these test functions, the right hand side of (|2.29l) is greater than or equal to 

5>|7 A|V$f +1 | 2 ^-E / ||S?|| L o 0(ni) |9 i $f+ 1 ||$f +1 |^ 
1=1 l Jn ' 1=1 Jn i 

+5/ i$n 2 ^+x: / *u'i*r i i a d* 

4jfi ' j'eJ, ,/r i.«' 
+ jf |^°° [(I + c) ef +1 - F(u1 +1 ) + F(u)] exp(-crf)di} gi tf +1 dx 

> X>{/ A|V$f +1 | 2 dx-f: / I^IUoo^lflj^H^ldx 

1=1 l/ n < l= l J ^ 



1$ 



Hi 



fe+1 12 



(2.30) 



> 



(| + c ) e ? +1 - + ^(«)] exp(-at)dt 



9m 



k+1 dx 



-,fc+l|2 





r r 










[(?- 






Jo 







> 



> 



g^/j,v W +§/„i*? 

p r /*oo "I r /»oo ^ 

+ J J (- + c - C'J ef +1 exp(-at)dt y ef +1 exp(-ai)dt g 2 do:j 



S>{/ 51 

1=1 lJf2 < 2 



where 



C" = C if J™ ef +1 exp(-at)di > 0, 



[ C' = -C\i J °° ef +1 exp(-at)dt < 0, 
and notice that a is large enough. 

Similarly, we can estimate the left hand side of (|2.29[) . which is in fact less than or equal 
to 



Wf A|V<&f||V^| C te+ / 2a\$t\\tf\dx 
l=1 Un t J tii 



J 

< 2 M i[ A (ll V *flli»(fl,) + ll V ^ 



fe||2 



£||$*M2 

2 



a n frn2 

2 1 '" ' 



+ E IIwmIU-ct,,,,) (11^11^(00 + iirfi 



(2.31) 



i 1 a J 1 



where Mi is a positive constant depending only on {^z}ie{i,j} and the coefficients of 
(|2.3[) . Since a can be chosen such that a > (maxj e {i,j} ||-B'|| ioo ,Q s) 2 , there exists M2 
positive, depending only on \S^i\i^{\,i} an d the coefficients of (|2.3p such that the right 
hand side of (|2.31[) is less than 



£m 2 [jT Qiv$f + 1 2 + ^iv$f +i i 2 + f i$r i 2 ) ^ 



(2.32) 



z=i 
Define 



< ?||V$f|| 2 2(ni) + -||$f|| 2 



A 



L 2 ^,) T 2ll V$ ! llL 2 (n,) 



fe+l||2 



1$ 



k+1 |2 



7 = 1 



? f lll 2 (a,) 
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lli^) I , 



(2.33) 



then from (pOUj) . (j2~3Tj) and (pT32j) . 

(/3 - M 2 )£ fe+1 < M 2 E k , 



(2.34) 



where /3 = min{/3i, . . . , /3/}. 

Since M2 depends only on {tli}^^ jy and the coefficients of (12.31) . (3 can be chosen large 
enough, such that 

M 2 

M 3 := 



/3-M 2 



< f . 



We obtain 



< M^Quv^ill^ + fll^lli^))- 



That implies 



||$ fe H 2 



I ML 2 (f2,) 



^ M 3 E f^l|V*?||i a( n,) + ll^ll^cn,)) • ( 2 - 35 ) 



Notice that (|2 . 35[) still holds if M3 and A are fixed, and a is replaced by all y which is 
larger than a. This observation leads to 



< MS 



Y,f Q (j[ ef exp(-yt)dtg^J dx 

4 \ 7 /• / z" 00 \ 2 

(j Q |Ve?|exp(-yt)tftJ ,g 2 dx 

+ -jllj n e?exp(-i/i)d^ |V. 9; | 2 ^ 

+ E / (/ e?exp(~yi)dt) gfds . 
j =1 \./o / 



(2.36) 



Let a' be a constant larger than or equal to a, we obtain from the previous inequality 
that 



< MS 



I r pCx'-\-X / roo \ 2 

Y. / ef exp(-yt)dt) gfdydx 

l=1 J Sit J a' \J0 / 

1 /■ /-a'+l / <.oo \ 2 

J2J n J, j[J \Ve°\exp(~yt)dtj gfdydx 



(2.37) 



7 r ra'+l ^ / roc \ 2 
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+ H / e?e>q>(-yt)di) sfdydx 

, =1 JO; ia' WO / 



Using the fact that u° belongs to C^°(il x (0, oo)), we can infer that the right hand 
side of (|2.37p is bounded by a constant M3 M±(a). Since gi is greater than 1, then 

Yl J J yj ef exp(— yi)dtj dydx < M^M 4 (a). (2.38) 

H2I55D infers 




(2.39) 



that concludes the proof. 

□ 

3. Convergence for Semilinear Elliptic Equations 

Consider the semilinear elliptic equation 

_ Xa7=i dj(a l . J d l u) + Yh=i b A u + cu = in 

(3.1) 

u = g, on <9f2, 

where fHs a bounded and smooth domain in R™ . We impose on the coefficients of (|3.1j) 
Conditions (Al), (A2) in the previous section and the following condition 
(A3') There exists C > 0, such that C < c(x) on H and for all x in E": 

\F(x,z) ~ F{x,z')\ < C\z- z'\,Vz,z' e E. 

With Conditions (Al), (A2) and (A3'), Equation (|3.1[) has a unique solution u in 
W l - 2 (U) nL°°(n) (see O, Hi). 

We impose the same way of dividing the domain VL and the same notations as in the 
previous section. 

The Schwarz algorithm at the iterate #k in the l-th domain, denoted by uf, is then 
defined by 

{ ~Yh n=i 9 j( a i,j d i u i) + E"=i kdiUi + cu\ = F(x,uf), in tt h 

< (3.2) 

{ = <b m >^~\ on i>,vz' e j,, 

where 05/^' are either Dirichlet or Robin transmission operators. 
Each iterate also inherits the boundary conditions of u: 

Ui = g on dfti D dft. 

A bounded initial guess u° in C°°(0 x (0, 00)) is also provided 

® M /«i=u onr, |i; ,vi'e Ji. 
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The difference ef between uf and u is a solution of 

- Tnj=l dj( a i,j d i e l) + b i d i e l + ce l = F ( X i u l) ~ F ( X i u )> in 

(3.3) 

58,,1/ef = B M /e?r\ on r,, r ,VZ' € Jj. 
Moreover, 

ef = on <9ft;. 

Classical Schwarz Methods 
By induction, each sub-problem (|3.2p in each iteration has a unique solution in W 1:2 (ft)D 
L°°(f2) for the Dirichlet transmission condition. The algorithm is well-posed. 
Consider (|3.3j) and let 5 be a bounded and strictly positive function in C 2 (R",R). Define 
the following function 

*?(*) := (ef (z)) 2 5 (z)- 
Let q be &^ + Y^=x % a i 1 jdj99~ 1 'j then belongs to L°°(Q;), and define 

n n 

£^($) = - £ a j -(o ili fi| e $)+ £cift$. (3.4) 

Lemma 3.1. Choose g to be g" 1 where g is a solution in C 2 (Q) nC(O) of the following 
equation 

j E"i=l d j( a i,o d iS) + 2(C+ ||c||oo)£f - Y!i=\ h i d i9 < 0, Wlfy, 

[ g is strictly positive and bounded on fl, 
then £;£i(<f>J £ ) < m £/ie distributional sense. 
Proof. Define the operator 

n 

£>WQ '■= — E dj{a<ijdi). 
A complicated but easy computation gives 

(n \ n 

- E *>,!<•,.,<>,< 1 i efg ~ J2 ^i,Aefd je fg - (3.6) 

n I n 

- E Aa hj d * e i e i d j9 + ( e f) 2 X 1 "•■i' , - :9 ' <>.,"■ j'^U) 
That implies 

n I n \ 

+ E 6 » + E toXjfygg- 1 < (ef ) 2 OT, (3.7) 

i=i V 3=1 / 
17 



(3.5) 



where 



m=\J2 ^(a iJ a lCg -i)) + 2 (c + || c || 0O )( 5 - 1 )-^M l ( 5 - 1 ) 

,t,j = l i=l 



9 2 - 



Therefore, the nonlinear equation (|3.3p has been transformed into the following linearized 
inequation of $ z fe : 

n n 

~ 9 M,A®l) + i)0i$? ^ °" ( 3 - 8 ) 

i,j=l i=l 

□ 

Theorem 3.1. Consider the Schwarz algorithm with Dirichlet transmission condition, 

lim max - u\\ L co(q ) = 0. 

Proof. Step 1: Construct some estimates of the errors {ef }. 
Consider (|3.8p and define 

A/ = esssup $f(x). (3.9) 

an, 

By the weak maximum principle, Theorem 8.1, [10(, $f(x) is bounded by M almost 
every where on f2;, that implies the following estimates, 



(ef) 2 g < max ^esssup (efY g J , a.e. in £1 ; , V/ e {1, . . . , /}. (3.10) 

Step 2: Convergence of the Algorithm. 
Denote 

E k = max (esssup (ef) 2 ^ . (3-H) 
From (|2.17[) . for all I' in Ji, and for a.e. x in Ti^ 

(e^(x)) 2 g{x) < max I esssup (e^r 1 {x) 2 g{x)) I , (3.12) 

l " eJ '' V r i',i" / 



or 



(ef(:r)) 2 < ^ mjix ^ SSU p(4- 1 (x) 2 3 (x))j (3.13) 

< g(x) max esssup (e| £ , _1 (x) 2 g(x) _1 ) 

Fix / in {1, . . . , /} and let / be a function in C 2 (H) such that 

18 



• / > on f2j/, W e Ji; 

• There exist two positive real numbers e small and M large such that for | V/(x) | < e, 
Yl .<),«,, JJ.fnr) >M; 

• / = on 0JV, Vr G J ; ; 

• 1 1 /I |oo — 1- (We can construct this function by constructing a function g which 
satisfies the fist three properties, and then take / = g/||g||oo)- 

Let p be a constant and put 

g = M - exp(pf) = exp(2p) - exp(p/), 

then, 



Y '>.,'■ 'i >.,'),>!) + 2(C + | |c| |oo)5 - ^2 b * d i9 



n 

+2(C + \\c\ |oo)(exp(2p) - exp(p/)) + J> eMpf)P^f 



< exp(-pf) 



-p 2 M + ep + 2(C + ||cj|co) 



/exp(2 i o) \ .., .. 

( v e^7)- 1 J +m f x||6l|Uep 



(3.14) 



< exp(-Af) -p 2 Af + ep + 2(C+||c|| 00 )exp(2p) + max 1 1 c>i 1 1 oo ep 

L i 

< o, 

when p is large enough and M > exp(2p). 
Inequality (I3.13[) then becomes 

(ef (x)) 2 < (M ~ exp(pf{x))) max I ess sup ((efr 1 ^)) 2 (M - l)" 1 ) ) .(3.15) 

'" eJ '' V r ''.'" J 

Since lies inside /(x) is strictly positive on !Ty/. Hence, there exists strictly 
less than 1, such that 

(ef (x)) 2 < Mi j max I esssup (ejT 1 ^)) 2 ] < Mi ,E fe_1 ,Vi' € J;, for a.e at in T ; /<.3.16) 

Similarly as in (13.26[) . let / be a function in C 2 (Vl) such that 

• / = on d$li\ 

• / > on on fi;; 

• There exist two positive constants e small enough and M large enough such that 
for \Vfi(x)\ < e, we have that d ji a i,j d if)i x ) > M ' 
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• ll/Hoo = l. 

and let p be a constant large enough. Setting 

g = Mo-exp(pf), 



(efOr)) 2 < s(sW(sup(e?(*)) 2 (Mo-l)) 
1 eJl \ r i,i> J 

< max I sup (ef(x)) ] , for a.e x in Qi 



(3.17) 



Combining (|3.16p and (I3.17[) . we can deduce that there exists M 2 strictly less than 
1 satisfying 

E k < M 2 E k ~\ (3.18) 

that leads to 



lim E K = 0. 

k— too 



□ 



3.2. Optimized Schwarz Methods 

The optimized Schwarz waveform relaxation algorithms are defined by replacing the 
Dirichlet transmission operators by Robin ones 

n 

*Bi,i>v = ^2 a i,j9ivni t i' j +pi,i'V, 

where riu'.j is the j-th component of the outward unit normal vector of T;.//; pij/ is 
positive and belongs to L^iTi^). By induction argument, each subproblem (|2.2I) in 
each iteration has a unique solution in H l (Vl) and the algorithm is well-posed. 

In general, optimized Schwarz algorithms do not always converge when applied to 
elliptic equations, as shown in the following example. 

Example 3.1. Consider the elliptic problem on the domain ft — (0, L) 

{ u" - 3u' - 4u = /, m ft, 

(3.19) 

{ u(0) = u{L) = 0, 

where f belongs to C°°([0,L]). Divide ft into two subdomains fti — (0,L 2 ) and ft 2 
/. i . /. : , with < L\ < L 2 < L, and consider the domain decomposition algorithm 

J (^ +1 )"-3K +1 )'-4^ +1 -/, in(0,L a ), 

1 ^ +1 (0) = and {u k+1 )'{L 2 )+pu k+1 {L 2 ) = (u k )> (L 2 ) + pu k (L 2 ), 
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r (^ +1 )"-3(^ +1 )'-4u 2 fc+1 = /, in{L x ,L), 

\ u k +\L) = and - «7^ +1 (ii) - - <ptf(Li), 

where p, q are positive numbers. 

The errors e\ and e 2 from the above equations can be obtained 



A k+1 (exp(4x) - exp(-x)), 



k+l 



B k+ i(exv(A(x - L)) - exp(-(a: - L))), 



wht 



T-2 



Set 



the 



A 



A'+l 



Bk+i 



4exp(4(L 2 - £)) + exp(-(£ 2 - L)) + p(exp(4(£ 2 - L)) - exp(-(£ 2 - £))) 
4exp(4L 2 ) + exp(-L 2 ) +p(exp(4L 2 ) - exp(-L 2 )) 
4exp(4Li) + exp(-Li) - g(exp(4Zi) - exp(-Li)) 
: 4exp(4(L 1 - L)) + cxp(-(Li - L)) - g(exp(4(L 1 - L)) - cxp(-(Li - L))) ' 



4exp(5L 2 ) + exp(5L) +p(exp(5L 2 ) — exp(5L)) 



4exp(5L 2 ) + 1 +p(exp(5L 2 ) - 1) 
4exp(5Li) + 1 - g(exp(5Li) - 1) 



(3.20) 



4exp(5Li) + exp(5L) — q(exp(5Li) — exp(5L)) 

The algorithm converges if and only if t is smaller than 1 . 
For p — 1 and q large, 

exp(5Li) — 1 



exp(5Li) + exp(5L) 



(3.21) 



Since the right hand side of <\S. 21\ is greater than 1 for L\ > L/h\2 and L large, the 
algorithm does not converge for p = 1 and q large. 
Note that 

4exp(5L 2 ) + exp(5L) + p(exp(5L 2 ) — exp(5L)) 



n 



4exp(5L 2 ) + 1 +p(exp(5L 2 ) - 1) 
can be made larger than 1 since 

4exp(5L 2 ) + exp(5£) 



4exp(5L 2 ) + 1 
is larger than 1; 

and the second term can be made larger than 1 since 



T~2 



exp(5£i) — 1 



cxp(5Li) + exp(5i) 



is larger than 1 for L\ > L/ln2 and L large. 
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Remark 3.1. In the above example, the Schwarz algorithms converge if p and q are 
large. This observation leads to Theorem \3.2\ below. Moreover, the convergence factor is 
still small if p, q are chosen to be negative. At least, the behavior of t is quite the same 
when q tend to +00 and —00. This observation is different from what was seen from the 
convergence rate in the case where the subdomains are two half lines in 0/ and @y. 

Introduce the modified Robin transmission operator, based on (T 2 ) 

n 

where p is a positive parameter. 

Theorem 3.2. Consider Schwarz Algorithms with Robin transmission conditions, if we 
replace 05;,;/ by 58?;, , then there exists po such that for p > pq, the algorithms converge 
in the following sense 

lim max \\uj — u||r,2m,i = 0. 

k^oo le{l,...,I] v ' 

Remark 3.2. Consider again Example 3.1, and a Schwarz algorithm which diverges with 
p = and q = qo, there exists L\, L2, L such that the algorithm diverges 

( {u k )"-Z{u k )>-Au k = f, m (0,L 2 ), 
1 u k (0) = and (u k Y(L 2 ) = (u*- 1 )'^), 

( (u k )" - 3(u k Y - Au k = /, m (L U L), 

\ u k (L) = and (u|)'(ii) - «fo«£(Li) - - qout^L,), 

where qo is a large constant. 

Let P be a function in C 1 (R), and put w — u exp(P), Equation (j 3. 19\i can be transformed 
into 

( w" - (3 + 2P')w' + (-4 - 3P' + {P'f - P")w = f, in (0, L), 
\ w(0) = w{L) = 0. 
The Schwarz algorithm then becomes 

( (w k )" - (3 + 2P')«)' + (-4 - 3P' + {P' f - P")w k = f, m (0, L 2 ), 
\ w k (0) = and ((w k Y - P'w k ){L 2 ) = ((w^ 1 )' - P'w^ 1 )^), 

( {w$)" - (3 + 2P / )(w k Y + (-4 - 3P' + (P') 2 - P")w k = /, in {L X ,L), 

\ w k {L) = and (w 2 fc )'(Li) - (P' + q )w k (L 1 ) = (K'- 1 )'^) - (P' + q )w k (L 1 ), 

We can deduce that given a pair of numbers (p, q) , we can find a class of functions P 
such that —P'(L,2) = p and —P'(Lx) = q + qo, and the Schwarz algorithm with the associ- 
ated equation t\3.22]i and this Robin transmission condition does not converge. However, 
Theorem \3.1\ announces that we can make the algorithms converge, even if they do not 
converge initially, by increasing the parameter p. 
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(3.22) 



Proof. Step 1: Linearize the equation (|2.3p . 

Consider the equation (|2 .3[) and let <?/ be a strictly positive bounded function in 
C 2 (£1/,R). Define the following function 

$?(aO :=ef(x) 5/ (x). 

A complicated but easy computation gives 

= - £ flif^j + X; + £ Oij U$f^f + fl,^ (3.23) 

i,j=i i=i i,i=i \ 9i 9i / 

+ ( £ - e 3^ - e + (« - Dtf )-) f . 

\i,j=i 5 ' i,j=i W > i=i 9 i J 

where 

f F(a;) = if uf(x) = u(x),x £ fl, 

I ^) = F( ^gi::& (g)) if«^(»)^u(»),»6n, 

F is then bounded as F is Lipschitz. 

Similar as in (|3.8p and in Step 2 of the proof of Theorem 13.11 we rewrite the last term 
on the right hand side of (|3.23[) into the following form 

n n \ 

- £ 9 i (a i , i 5 i (( 9 f)- 1 )) + (c-F)(( 5 f)- 1 )+^6 i 5 i (( ff f)- 1 ) Uf, (3.24) 

i,i=l i=l J 

and use the same argument as in Q3.5p : choose g\ to be g; -1 where gi is a solution in 
C 2 (f2) fi C(Q) of the following equation 

- Y%j=i di(ai,j9i§i) - Kgi + Yh=i b i 9 i9l > 0, in Qj, 

_ (3.25) 
g is strictly positive and bounded on f2, 

where K is a. large enough constant. 

Since pi i> is strictly positive for all I in {1, . . . , 1} and /' in J;, there exist functions 
fl, I £ {1, . .'. , /}, in C 2 (f2) such that 

• S"i=i "...m.r.j^-fi =Pl,V on r /,;', V/' G J;. 

• S"i=i "'.."r.i.j'lj: =Pl>,l on IY/, V/': Z G Jv . 

• There exist two positive constants e small enough and M large enough such that 
for |V/i(as)| < e, £!\ =1 ^-(a^/)^) > M. 

• fi = 0, on IYp, VZ' G J;; and /// = a u on \/V : I £ Jy. 

• 1 1/| |oo = I- (We can construct this function by constructing a function g which 
satisfies the fist properties, and then take / = ff/Hslloo)- 
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Similar as in p.26[) . let p be a constant large enough and put g = M3 — exp(— pf), where 
M3 is a positive constant, 

n n 

^2 dj( a i,jdi(j) + 2(C + ||cl|oo)g - y^bjdjg 

n n 

= - E P 2 exp(pf)d j {a iJ d t f) - ^ <H,jPexp{pf)9i,jf 

n 

+2(C+\\c\\ 00 )(M 3 -exp(pf)) + Y / b l exp(pf)pd t f 

i=l 

n 

-Xp 2 M - 2J cti,jdi,jfp 

+2(C+|| c |U) M3 e ; p e ( y ) + ex P (p/) 
< 0, 

when p is large enough. 

Denote the right hand side of the equation (|3.23[) by £/($ z fe ), then it can be rewritten 
in the following form 

n n 
— 1 i—1 

where B\ and C l are functions in L°°(M. n ), C l is bounded from below by K + c + F. p 
can be chosen such that there exists a large enough, 2a > Ci > a. 

Now, consider the Robin transmission condition on the boundary Fiji 

n 

®i(*j fc ) = E ».J>M»<r.j ( 3 - 27 ) 

E a hj n U',3 d i e l } 91 + E a i,3 n l,l',j®i9l e l 
\i,j=l j »,i=i 

\ 

E a i,j n l',l,jdi e l +Pl,l' e l 9l 

v»,i=i / 



E a i,j n l',l,3 d i e l' 1 +Pl,l' e l< 91 
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We can choose /; such that 



- Pi, on Tu,,W e Ji, 



fji 

where /3i is a constant greater than 1. 

From the previous calculation on Q3;(<I>f) and £;($j), $>f is in fact a solution of the 
following equation 

f = 0, inO ; x (0,oo), 

i (3.28) 
1 /3i«8i,,/($f) = Q3ji'($f, _1 ) on r,,,, x (0,oo),W e Jj. 



Step 2: The Proof of Convergence. 

Denote by Q; to be the open set f2i\U;/ e j i f2//. For each Z in {1, . . . , /}, let <pf to be 
a function in 1 (O; ) and <pf +1 to be a function in such that <^f +1 = <pf, on IY;' 

for all /' in J;. Now, using <ycf +1 and c/?f as test functions for all subdomains, we obtain 



I [ ~ n » n » 

- E / pi',i®i<pi da \ 

i=i ij=i "' n i i=i 

+ / C'$f+Vf +1 ^+ X) / PM' $ f + Vf +1 ^|- (3-29) 

JOl 2'£jj'' r M' ) 

In the above equality, choose <^f +1 to be then there exists (^f such that = 

£^^, +1 on T;.;' for all i' in J;; and 

With these test functions, the right hand side of (|3 . 29[) is greater than or equal to 

5>{/ A|V$? +1 | 2 ^-£ / ||^|Uoo (n() |a i $f +1 ||$f +1 |da ; 
;=i i=i 

+af |$f +1 | 2 ^+X: / Pr,^f +1 | 2 ^) 



> X> 
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-a / |$f +1 | 2 



i=1 L/n, ^ * J a, 



with a being large enough. 

Similarly, we estimate the left hand side of (|3. 291) . which is in fact bounded by 



Wf A|V*f||V^|dx+ / 2a\^\\^\dx 
l=1 UOj Jiii 

n p 



illL°°(f2,) 



I'eJi r i'.i > 



< 



a||$ fc " 2 



^M 4 [A(||V$f|| 2 2(f2() + ||V^|| 2 

z=i 

+ (llVSfll^ + ( max } |N|| ioo(S() ) 2 ||^|| 2 i2(fti) 
£ llw,'IU~(r,/,,) (ll $ 



L»(n,) 



L2(n ; ) 



(3.30) 



where M4 is a positive constant which depends only on and the coefficients 

of (|3.3[) . Since a can be chosen such that a > (maxigj^ . ^) II^IIl°°(si ; )) 2 ' there exists 
M5 positive, depending only on {£!/};£{! n and the coefficients of (|3.3p such that the 
right hand side of (|3.30p is less than 



Em 5 [/ (W| a + f|* 



Z=l 



'• 2 |_ ^|v$f +1 | 2 + -|$f +1 | 2 



dx 



(3.31) 



z=i 



Define 



A 



llL 2 (0,) 1" ^ll^i llL 2 (f2() I J 



fc||2 



(3.32) 



(3.33) 



then from (|3T5D]) . ([331]) and 

(/3-M 5 )£ fe+ i < M 5 £ fe , 

j8 = min{/?i, . . . , Pi}. Since M5 depends only on {£!;};ei.j and the coefficients of f|2 . 3[) 
/3 can be chosen large enough, such that p M l Ir < 1, then 



Mr 



fe-1 



Si, 



which means £Jfc tends to as fc tends to infinity. This concludes the proof. 
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(3.34) 
□ 



4. Conclusions 

We have introduced a new class of techniques to study the convergene of Schwarz 
methods. In particular, classical Schwarz methods are proved to converge when being 
applied to both parabolic and elliptic equations. On the contrary, Schwarz methods with 
Robin transmission conditions only converge when we use them for parabolic equations, 
though they were proved to converge faster than classical ones in previous studies. For 
elliptic equations, we have given a counter example, where we can see that optimized 
Schwarz methods do not converge; and for each optimized Schwarz algorithm, there ex- 
ists a class of elliptic equations which is not stable with this algorithm. A new way of 
stabilizing the algorithms has then been proposed. 
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